Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



> 



Star-forming protoclusters associated with methanol masers 

V. Mimeri'2, M.G. Burton^, T. Hill^, M.R. Pestalozzi^, C.R. Purcell^, G. Garay"^, A. Walsh^ and 

S. Longmore^ 

■ ^ Service d'Astrophysique, DAPNIA/DSM/CEA CE de Saclay, 91191 Gif-sur-Yvette, France 
, e-mail: Vincent.Minier@cea.fr 

Q ■ ^ School of Physics, University of New South Wales, Sydney 2052, NSW, Australia 

^Nj ' ^ Onsala Space Observatory, S-439 92 Onsala, Sweden 

►-^ , * Departamento de Astronomi'a, Universidad de Chile, Casilla 36-D, Santiago, Chile 

Q ■ ^ Center for Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA 

' Received / Accepted 

<N : 

' Abstract. We present a multiwavelength study of five methanol maser sites which are not directly associated 

' with a strong (> 100 mjy) radio continuum source; G 31.28+0.06, G 59.78+0.06, G 173.49+2.42 (S231, S233IR), 
G 188.95+0.89 (S252, AFGL5180) and G 192.60-0.05 (S255IR). These radio-quiet methanol maser sites are often 

' interpreted as precursors of ultra-compact Hll regions or massive protostar sites. In this work, the environment 

, of methanol masers is probed from mid-IR to millimetre wavelengths at angular resolutions of 8" — 34". Spectral 

^— ( ' energy distribution (SED) diagrams for each site are presented, together with mass and luminosity estimates. 

T— I ' Each radio-quiet maser site is always associated with a massive (> 50 M0), deeply embedded (Av > 40 mag) and 

. very luminous (> 10* Lq) molecular clump, with LtotaiocMgaJ' ■ These physical properties characterise massive 

' star- forming clumps in earlier evolutionary phases than Hll regions. In addition, colder gas clumps seen only at 

fH I mm- wavelengths are also found near the methanol maser sites. These colder clumps may represent an even earlier 
phase of massive star formation. These results suggest an evolutionary sequence for massive star formation from 

I ' a cold clump, seen only at mm wavelengths, evolving to a hot molecular core with a two-component SED with 

O , peaks at far-IR and mid-IR wavelengths, to an (ultra-compact) Hll region. Alternatively, the cold clumps might 

■ be clusters of low-mass YSOs, in formation near the massive star-forming clusters. Finally, the values of the dust 
^ ' grain emissivity index (/3) range between 1.6 and 1.9. 

^ ' Key words, masers — stars: formation — circumstellar matter 
•1— I . 

^ 1. Introduction in fact generally not co-spatial with strong radio sources 

(e.g. UC Hii regions), but frequently tend to be isolated 
Class II methanol masers were originally detected toward ^^^^ ^^^^^^ Out of ~ 250 methanol maser sites stud- 
regions of massive star formation (MSF) because many -^^^ 75% ^-^^^^^^ associated (within 0.1 pc) with 
searches at 6.7 and 12.2 GHz had focused on observmg Hii regions or more extended radio sources (Caswell 
regions exhibiting typical MSF signatures, such as strong ^^^^^ pj^jUip^ y^-^-^^ 
radio continuum sources (i.e. Hii regions), IRAS sources y^^y ^^^^ 4 ^jy) ^^^-^ continuum emis- 
with IR colours of ultra-compact H II (UC H ll) re gions ^j^^ however, been detected toward isolated methanol 
and OH masers (e g. Batrla et al.[M7| Schutte et al.|MI maser sites (e.g. van der Walt et al.m^- An extreme case 
CasweU et al.im. At low angular resolution (- 1'- 10'), iraS 20126+4104 where the methanol maser in the hot 
many 6.7 and 12.2-GHz methanol masers appeared to be ^^oiecular core is associated with a 0.1-mJy radio source 
associated with H 11 regions. g g qjj^ (Hofner et al. lMTHI Minier et al. lWlTll . This 
In contrast, high angular resolution (~ 0.01" - 1") ob- might be the first sign that some isolated methanol masers 
servations have demonstrated that methanol masers are are associated with precursors of UC Hii regions, i.e. pro- 

tostellar objects on the way to forming a star or a multiple 

Send offprint requests to: Vincent Minier g^^r system. Methanol maser sites are therefore generally 

' In this paper, a UC H 11 region is defined as a region ionised ^adio-qmet in terms of continuum emission (< 1 mJy at 5 

by high-mass stars, with a diameter ~ 0.1 pc and an electron , . j_ ■ 1 u i- i -i- t 

, . ,„4 i , . ^ r kpc, I.e. a typical observational sensitivity limit), 

density >10* cm"^ producing a 22-GHz flux of ~ 100 mJy at ^ ' •'^ J ) 

a distance of 5 kpc (Kurtz et al. 
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Complementary spectral line observations confirm this 
view. Preliminary results from a molecular line survey 
with the Mopra millimetre telescope suggest that isolated 
methanol masers are possibly associated with hot molec- 
ular cores (HMCs) (Purcell et al. '2004!'), which are pre- 
sumed to be the sites where massive protostars evolve to 
form an ionising ZAMS star (Kurtz et al. 2000; Fontani et 
al. 2003). The intense and complex chemistry, the temper- 
ature (some 100 K) and the density (>10'' cm~'^) of these 
HMCs would provide suitable conditions for methanol 
masers to arise (see Minier et al. 2003 and references 
therein). 

Finally, recent ATCA observations by Minier et al. 
H2003|l might indicate an exclusive association between 
methanol masers and massive star-forming regions follow- 
ing a nil detection rate of masers toward regions of low- 
mass star formation. All these elements tend to demon- 
strate that methanol masers may trace an early stage of 
MSF. Whether methanol masers trace embedded massive 
protostars or suitable conditions nearby a site of MSF re- 
mains an open question. 

This paper is devoted to a multiwavelength study 
of five radio-quiet methanol maser sites and their 
environments, which are particularly good candidates 
for hosting a massive protostar (Minier et al. 120001 
TMrn -. G 31.28-1-0.06, G 59.78-1-0.06, G 173.49-1-2.42 (S231, 
S233IR), G 188.95-1-0.89 (S252, AFGL 5180) and G 192.60- 
0.05 (S255IR). Comparisons of images taken from optical 
to radio wavelengths are presented, accompanied by es- 
timates of the mass, luminosity and temperature of the 
MSF clumps in each region. We will demonstrate that the 
five radio-quiet maser sites coincide with massive and lu- 
minous star-forming clumps. Finally, it will be shown that 
methanol masers trace very young clusters of star forma- 
tion in the earliest stages of their evolution. 

Sect. 2 describes the observations and data archives 
used for this work. Derived physical quantities and di- 
agrams of spectral energy distributions are presented in 
Sect. 3, followed by a discussion in Sect. 4 and conclu- 
sions in Sect. 5. 

2. Observations, data archives and analysis 

2.1. Optical image archives 

An optical image of an 8' x 8' region was generated 
around each source using the Digitized Sky Survey (DSS) 
database (http://archive.eso.org/dss/dss). The Digitized 
Sky Survey comprises a set of all-sky photographic surveys 
conducted with the Palomar and UK Schmidt telescopes. 
The objective was to visualise any association with optical 
H II regions or main sequence OB stars after identification 
with the SIMBAD catalogues on CDS. 

2.2. Mid- and far- infrared image archives 

The mid-infrared (mid-IR) and far-infrared (far-IR) im- 
ages were obtained from the Midcourse Space Experiment 



(MSX) and from the InfraRed Astronomical Satellite 
(IRAS) data archives. The values for the mid-IR fluxes 
were estimated from the MSX images using the KARMA^ 
package, by estimating the flux in the source and also in 
the background. The background residual flux was mea- 
sured by taking the averaged flux in boxes at diverse posi- 
tions around the mm source. This method allows removal 
of the background noise as well as emission from nearby 
objects and extended regions. To avoid confusion with sil- 
icate emission/absorption at 8.7 /im (in the MSX A-band) 
and with the more extended PAH contribution at 7.7, 
8.6 and 11.3 /im (in the A- and C-bands), the emission 
contours in the E-band (21.3 /im) were used to estimate 
the extent of the thermal dust emission. The source flux 
then had the background subtracted, and the value was 
converted from W m^^ sr^^ to Jy. The image sensitivity 
varies from 1.9 to 29.8x10^*" W m^^ sr^^, with an angu- 
lar resolution of 18" (i.e. 0.18 pc at 2 kpc). The positions 
in the MSX point source catalogue are accurate to within 
1" - 5". The conversion factors were 6.84x10^, 2.74x10-*, 
3.08x10-* and 2.37x10"* Jy per W m'^ sr"* for the A- 
band (8.3 /im), C-band (12.1 /im), D-band (14.6 /tm) and 
E-band (21.3 /im), respectively (see Egan et al. 1 19991 for 
a MSX explanatory guide). These factors include a cor- 
rection turning the 6-arcsec square pixel into a Gaussian 
area. The far-IR fluxes were taken directly from the IRAS 
Point Source Catalogue (IRAS PSC). 

2.3. Submillimetre and millimetre continuum 
observations 

The submillimetre observations at 450 and 850 /im were 
conducted using the Sub-mm Common User Bolometer 
Array (SCUBA) on the James Clerk Maxwell Telescope 
(JCMT) in September 2002 for three of the five sources. 
G31.28-h0.06 was observed by Walsh et al. ipUSjl . The 
observations were made in jiggle map mode with a map 
size ~2-3 arcmin. The pixel size was 3 arcsec. Pointing and 
calibration were made on HL Tau. The average zenith at- 
mospheric optical depths were 1.8 and 0.3 at 450 and 850 
/tm, respectively. The volt-to-Jy-beam~* conversion fac- 
tors were 373.1 and 260.0. The image sensitivities were 
^ 0.3 and 0.01 Jy beam"* at 450 and 850 /tm, respec- 
tively. The FWHM hea,m sizes were estimated to be 8 and 
15 arcsec (i.e. 0.08-0.15 pc at 2 kpc) at 450 and 850 /tm, re- 
spectively. The data were reduced with the standard pro- 
cedure SURF (reduce-switch, flatfield, extinction, remsky 
and rebin) . The flux density values were estimated in each 
image with the KVIEW procedure described in Sect. 2.2. 

The 1.2-mm observations were made using the SEST 
IMaging Bolometer Array (SIMBA) on the Swedish-ESO 
Submilhmetre Telescope (SEST) in October 2001 and 
June 2002. SIMBA is a 37-channel hexagonal array in 
which the HPBW of a single element is about 24" (i.e. 
0.24 pc at 2 kpc) and the separation between elements on 
the sky is 44". The bandwidth in each channel is about 

^ See: |iittp://www.atnf.csiro.au/computing/software/karma/ 
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50 GHz. Spectral line emission present in the band may 
contribute to the total continuum flux values reported in 
this paper, though not more than 30%. Free-free emission 
is not expected to contribute to the (sub-)mm flux given 
the radio-quiet nature of the targets. Assuming optically 
thin free-free emission for frequencies > 8 GHz, 80 mJy 
at 8.6 GHz would convert to 57 mJy at 250 GHz, i.e. 
~ 1% of the 1.2-mm flux measured at the same position. 
Similarly, for an ionized wind model with (Panagia 
& Felli 11975(1 . the extrapolation at 250 GHz of the ra- 
dio fluxes would yield less than 11% of the observed flux 
densities. These observations were made in fast-mapping 
mode. Typical areas of 600" x 384" were imaged in or- 
der to detect nearby IRAS sources in the fleld, using a 
scan speed of 80" s~^. The resulting pixel size is 8 arcsec. 
Typically, total observation times were around 15 min- 
utes per source, giving a sensitivity 0.1 Jy beam"^. 
Zenith opacities ranged from 0.2 - 0.4 throughout the ob- 
serving runs in 2001 and 2002. The resulting data were 
reduced with the MOPSI mapping software package de- 
veloped by R. Zylka, using the "deconvolution" algo- 
rithm'' to remove the contribution of the electronics aris- 
ing from the fast-mapping observing mode, the "convert- 
ing" algorithm (Salter 11983(1 to convert the coordinates 
from rectangular to equatorial, and partly the N0D2 and 
GILDAS libraries. The calibration of our data was per- 
formed using Uranus. The multiplication factor between 
counts and Jy beam~^ was 138.8 mJy count"^ beam^^ 
for the October 2001 run (G 188.95-h0.89, G 192.60-0.05) 
and 66.9 mJy count^^ bearn"^ for the June 2002 run 
(G31.28-f0.06, G59.78-h0.06). The flux density values were 
estimated with the KVIEW procedure described in Sect. 
2.2. The 1.2-mm fluxes for G 173.49-^2.42 were taken from 
Beuther et al. l|M7^ . 

2.4. Molecular spectral line observations 

The molecular line observations were carried out in May 
1999 with the Onsala-20m millimetre telescope (OSO- 
20m) and in October 2002 with the Mopra millimetre 
telescope. CH3CN (92 GHz) and C^^O (109 GHz) lines 
were observed toward each methanol maser position with 
a single pointing. 80 MHz and 64 MHz bandwidths were 
used at Onsala and Mopra, respectively. The gain of the 
Onsala-20m telescope is ^ 20 Jy K^^ in the range 85- 
112 GHz while the gain of the Mopra telescope is about 
30 Jy K~^. CH3CN observations were performed with 
Mopra only. C^^O observations were carried out with 
Mopra toward G 192.60-0.05 and with the OSO-20m to- 
ward the other four sources. The pointing errors were es- 
timated to be less than 5 arcsec for the OSO-20m and less 
than 10 arcsec for Mopra. The data were processed with 
XS, the Onsala data reduction software. 



^ Developed within the SIMBA collaboration. 



2.5. Radio continuum data archives 

Radio continuum information was obtained from Aladin 
and SIMBAD on CDS within a 4' x 4' region around each 
methanol maser site. Speciflc radio survey catalogues were 
also searched using VizieR. These included the 1.4-GHz 
NRAO VLA sky survey (NVSS, Condon et al. lTMSjl . the 
5-GHz VLA survey of the Galactic plane (Becker et al. 
11994(1 . the survey of radio H 11 regions in the Northern Sky 
through recombination line emission 

(LockmanHnnni), the 

GB6 catalogue of radio sources at 4.85 GHz (Gregory et 
al. ll996(l . the UC H 11 region survey (Kurtz et al. 1994) and 
an ATCA survey at 8.6 GHz (Walsh et al.^MB- Source 
positions from SneU & Bally ijTOSSIl for G 192.60-0.05 and 
from Sridharan et al. for G59.78-f0.06 were also 

plotted. 

3. Results 

3.1. Derivation of physical quantities 

The goal of this paper is to derive the physical characteris- 
tics of the methanol maser sites. These physical quantities 
include the luminosity {L total), the gas mass (Mgas), the 
gas column density ( A^H2 ) 1 the gas density (riHa ) , the dust 
(Tcoid and That) and gas temperatures {Trot), the virial 
mass (Mmr) and the power law index of the dust opti- 
cal depth (/?). The values of these physical quantities are 
derived from the dust continuum emission and from the 
molecular line emission for each clump in the five MSF 
regions studied in this paper. We do not intend to model 
density and temperature structures in detail. This will be 
undertaken with higher angular resolution instruments in 
the future. 

3.1.1. Spectral Energy Distribution of dust emission 

The Spectral Energy Distribution (SED) of each millime- 
tre source with continuum data points from 8 /im to 1.2 
mm is fitted with two grey-body functions modified with 
a non-constant emissivity function, which model emission 
from a hot (> 100 K) dust core (diameter, dhot) embedded 
in a larger cold (< 30 K) or warm (30 < T < 100 K) enve- 
lope (diameter, dcoid)- The dust optical depth in the emis- 
sivity function is expressed as a power law Td — tq{v/vq)^ 
where v is the frequency and (3 is the dust grain emis- 
sivity index. The hot dust core is characterised by a dust 
temperature {Thotj and a solid angle {Vlhot) which are two 
free parameters in the first grey-body model. The cold 
(or warm) envelope is characterised by Tcoid and Vlcoid in 
the second grey-body function. Two additional parame- 
ters, the turn-over frequency between optically thin and 
thick regimes [v—Vf) where t=tq=1) and /3, are also esti- 
mated. In total six free parameters are therefore estimated 
using a minimum x^-Levenberg-Marquardt bestfit of the 
functions to the data. 

This model is a good approximation if each methanol 
maser site harbours a single young stellar object (YSO). 
In the case of a multiple YSO system or a YSO cluster. 
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our model will reproduce emission from both the indi- 
vidual YSOs and the cluster envelope. However, only a 
few massive YSOs are expected to be present within each 
methanol maser site due to the relatively small region ob- 
served. The SED will be dominated by emission from the 
most luminous objects present. In summary, our model 
will give the SED of a single YSO at best and the total 
SED of a YSO cluster at least. 

The total luminosity {L total) is derived by integrating 
the SED over the frequency /wavelength range. The gas 
mass is estimated using the following equation, adopting 
the Rayleigh-Jeans approximation to the Planck function, 
and assuming optically thin emission at 1.2 mm: 



Mgas = 5 X 10 



2kTcoidRdKd 



in Mr^ 



(1) 



where i^i.2mm is the 1.2-mm continuum intergrated flux 
in W m~^ Hz~^, D the distance to the source in m, A the 
wavelength in m, k the Boltzmann's constant in J K~^, 
Tcoid the dust temperature in K, Rd the dust-to-gas mass 
ratio and Kd = kq{v/vq)^ the mass coefficient absorp- 
tion (ocTrf) in IT? kg^^. The Rayleigh-Jeans approximation 
holds as long as the temperature is greater than 40 K at 
1.2 mm or 250 GHz. For temperatures less than 30 K, the 
use of the full Planck function is recommended. The values 
of Kd and Rd may vary with physical conditions. Values 
of 0.1 m^ kg^^ at A = 1.2 mm and 1% are adopted for 
Kd and Rd, respectively, in the present work (Ossenkopf & 
Kenning 1 1994|l . This gives values of Mgas similar or down 
to 4 times less than when using the mass estimate method 
from Hildebrand H1983|l . 

A'hs and riH2 are calculated by dividing Mgas by 
7r(0sest-D/2)^ and ATT{9sestD/2)^ /3, respectively, where 
9sest ~ 24" is the angular resolution of SEST at 1.2 mm. 
The values of the column density and density in table 5 
are therefore beam averaged values. and the H2 mass 
density are also estimated using dcoid, which is the fitted 
diameter of the protostellar envelope (tables 5 and 6) . 



Using the equation derived in Nummelin (1998) and 
the above assumptions, the rotation diagram is a plot of 
the left hand side quantity in Eq. 2 versus E^/k. 



\hc^Auigu 



Thdv I = / 



N 

W) 



En 



hTr, 



(2) 



where Vui is the frequency of the rotational transition, h 
is Planck's constant, c is the speed of light, A^i is the 
absorption coefficient, gu is the degeneracy of the upper 
energy level, N is the total CH3CN column density and 
Q{T) is the partition function. If the hypotheses were cor- 
rect, the data points would fall along a straight line (see 
Fig. 1) whose inverse of the negative slope is the rotational 
temperature (Trot)- K all transitions are thermalised, the 
derived temperature is an estimate of the kinematic tem- 
perature of the molecular gas. Note that Ti, is calculated 
from T|, the antenna temperature as follows: 



rimbVbf 



(3) 



where rjmb is the main-beam efficiency factor and rjbf is 
the beam- filling factor. In this work, we assume rjbf = 1 
(beam-averaged) i.e. T^b = Tb and -qmb — 0.55 for the 
OSO-20m and ri„ib = 0.7 for Mopra. 

The C-'^^O line is a very useful tool for estimating the 
H2 column density. In the optically thin regime the C^^O 
column density is given by Bourke et al. H1997fl : 



Nqiso = 2.42x10^ 



0.88 



1 - exp(-5.27/Te,) 



J(Te,) - JiTbg) 



T„.b(C^«0)di; in cm-^ (4) 



where J(T) = hv / k{exp(hi' / kT) — 1), T^^ the excitation 
temperature, Tbg the background temperature and h is 
Planck's constant. 

Once Nqisq is known, N^^ is derived using the column 
density relationship in Frerking et al. I|1982B: 



3.1.2. Molecular line analysis 

The molecular lines in the CH3CN and C^^O spectra are 
fitted with Gaussian functions. The Gaussian parameters 
are given in tables 3 and 4. CH3CN(5if — Ak) spectral 
lines corresponding to different if-components are used to 
derive an estimate of the gas temperature using the rota- 
tional diagram method and assuming that CH3CN emis- 
sion is optically thin. The rotation method is based on 
three assumptions. (?) The spectral line is optically thin. 
(a) The brightness temperature (Tbg) of the background 
is negligible with respect to the brightness temperature 
{Tb) of the source, {in) The population of all observed en- 
ergy levels is described by a single Boltzmann temperature 

Trot • 



iVH2 = 5.9 X 10'^A^ci«o + 1-3 X 10^^ in cm"^ 



(5) 



LTE is assumed and T^x — Trot- In this way, the rota- 
tional temperatures deduced from the CH3CN data is used 
to derive the C^^O column density. The deduced H2 col- 
umn densities are representative for the CH3CN emitting 
regions if their angular sizes are similar to that of C-'^^O 
emitting regions. These assumptions will be discussed fur- 
ther in Sect. 4. The total H2 mass is derived by multiply- 
ing 7Vh2 by n{6osoD/2)'^, where 9oso ~ 34" is the angular 
resolution of OSO-20m at 109.8 GHz. The density is es- 
timated by dividing the total mass by a source volume 
V = 4:/3tt{9osoD/2)^. Mgas and na^ are beam averaged 
values, i.e. a lower limit of the actual values if the C^*0 
cloud has an angular size lower than 9oso- 



V. Minier et al.: Star- forming protoclusters associated with methanol masers 



5 



G31.28-MM1 



a 20 40 



100 120 140 



G173.49-MM1 



10 20 



30 40 50 60 70 

E Ik (K) 



G188.95-MM1 



10 20 



30 40 50 60 70 

E Ik (K) 



G192.60-MM2 



10 20 

Fig. 1. Rotat 
the left-hand 



ional diagrams where lii{Nu/gu) represents 
side of Eq. 2. 



3.1.3. Virial mass 

The virial estimate of the molecular mass follows the ap- 
proach by MacLaren et al. H1988|l with density profile in 
R~^, where R is the radius. The virial mass is given by: 



M„ 



= 126At/i? in Mr. 



(6) 



where a is the full 3-dimensional velocity dispersion, R is 
the radius of the cloud in pc, Av is the FWHM of the 
C^^O line in km s^^ and G is the gravitational constant. 
In this work, Aw = Awcisq and R = 6osoD/2. A molecular 
clump is gravitationally bound as long as Mgas/Myir > 1. 
Due to uncertainties in the mass estimate, the relation 
Mgas/Myir > 1 will be discussed for each object. Note 
that the values of Myir in table 5 should be divided by 1.4 
(= 34"/24") before comparison with Mgas derived from 
the dust continuum observations. 



3.2. Presentation of the results 

For each methanol maser site, an overview of the region 
consisting of a (sub)mm contour map and symbols rep- 
resenting radio sources, IRAS and maser sources overlaid 
on an optical image is given. This is image a. in Figs. 2- 
6. A close-up image of the methanol maser site is also 
presented with the mid-IR and 450 fiui contours over- 
laid on the 850 /im image. This is image b. in Figs. 2- 
6. For the mm emission sources with enough data points 
from mid-IR to mm wavelengths (table 1), the Spectral 
Energy Distribution (SED) diagrams are presented with 
the physical parameters and derived physical quantities. 
IRAS fluxes were used only when an association with a 
mm source was clearly identified (table 2). These are plots 
c. in Figs. 2-6. Finally, CH3CN and C^^O spectra are pre- 
sented in Figs. 2-6 d. with corresponding data values in 
tables 3 and 4. Rotational diagrams are given in Fig. 1. 
Table 5 summarises the derived physical quantities for the 
five methanol maser sites. To derive these values, the OSO 
and SEST telescope beams are mainly used as physical an- 
gular sizes because C^*0 line emission and dust emission 
may not arise exactly from the same region. An estimate of 
the clump diameter in the SED model is given in Fig. 2-6. 

3.2.1. G 31.28+0.06 

Three mm sources are detected toward G 31.28-1-0.06 (Fig. 
2a). They are distributed along a dusty filament seen 
with SIMBA (Fig. 2a). The methanol maser site coexists 
with a H II region and an IRAS source in the mm source 
G31.28-f0.06-MMl (hereafter G 31. 28-MMl), but is sepa- 
rated from the radio source peak emission by 3.5 x lO'' AU 
at a distance of 5.6 kpc. The morphology of the H 11 region 
is irregular and clumpy, with a major axis of about 10". 
From the values given by Kurtz et al. I)1994|l . the peak 
brightness temperatures are only ~ 80 K at 3.6 cm and 
~ 55 K at 2 cm. This implies that this Hii region has 
a low value EM, and does not belong to the class of UC 
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Table 1. Continuum fluxes. Positions are the 1.2-mni peak emission positions. Mid-IR fluxes are from MSX images. 
Sub-mm fluxes from SCUBA images. Mm fluxes from SIMBA images. means that no emission was detected. Notes: 
*fluxes interpolated from SED diagrams; **MAMBO mm fluxes from Beuther et al. (12002,1 . References: ^ Walsh et al. 
ifTOMjl : ^Kurtz et al. ifTUMjl : ^Sridharan et al. PUU^ : "Snell & Bally 1(1^ . 
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Table 2. IRAS fluxes and association with mm sources. 



Source 


IRAS name 


12 


25 


60 


100 


log(F25/Fi2) 


log(Feo/Fi2) 


Association 






(^m) 


(^m) 


{lim) 


{lira) 








G 31.28-^0.06 


18456-0129 


4.85 


89.33 


1071.0 


3693.0 


1.26 


2.34 


G31.28-MM1 


G 59.784-0.06 


19410-F2336 


14.4 


108.8 


982.5 


1631.0 


0.98 


1.83 


G 59.78-MM2 


G 173.49-1-2.42 


05358-H3543 


5.6 


74.7 


722.3 


1310.0 


1.12 


2.11 


G 173.49-MM4 


G 188.95+0.89 


060584-2138 


13.9 


140.2 


955.5 


1666.0 


1.00 


1.83 


G 188.95-MMl 


G 192.60-0.05 


06099+1800 


107.2 


371.6 


3145.0 


5285.0 


0.54 


1.47 


G 192.60-MM2 



Hii regions as defined by Kurtz et al. H2000|l . The mun- 
ber of ionizing photons required to excite the H ii region 
is ~ 4.6x10'*'' s~*. This could be supphed by an 09.5 
ZAMS star, which has a luminosity of 3.8 x 10^ Lq. No 
OB association or optical Hii region is seen nearby to 
G 31.28+0.06 although W 43 is located 25' south. The na- 
ture of the NVSS radio source in N-E is unclear. G 31.28- 
MM2 and -MM3 are much weaker mm emitters. G 31.28- 
MMl and G31.28-MM2 also coincide with sub-mm and 
mid-IR sources within the positional accuracy (Fig. 2b). 
In figure 2c, a SED model is presented for the clumps 1 
and 2. The SED of G 31.28-MMl is consistent with that of 
~ 10^ L0 and - 10^ Mq clump. The SED of G 31.28-MMl 
is not well modelled between 8.3 and 12.1 /xm. This may be 
due to silicate absorption at 8.7 /im that is not taken into 
account in our modelling. CII3CN and C^^O lines (Fig. 
2d) were detected toward the methanol maser site at ve- 
locities falling in the methanol maser velocity range (106 
to 114 km s~*). The gas column density and temperature 
agree quite well with those derived from dust observations. 
The clump is gravitationally bound (see table 5). G 31.28- 
MM2 is less luminous (10^ Lq), less massive (340 Mq) 
and much cooler (16 K) than G 31.28-MMl. 



3.2.2. G 59.78+0.06 



G 59.78+0.06 is a typical isolated maser, with no bright 
UC Hii region associated with it (Fig. 3). The maser 
site, however, coincides with an IRAS source satisfying 
the Wood & Churchwell criteria for a UC Hii region, as 
well as a weak radio source of flux ~ 1 mJy (Sridharan et 
al. l2()02|l . The 1.2-mm continuum map reveals a "peanut" 
shape with two clumps, G 59.78-MMl and -MM2. G 59.78- 
MMl is the coolest clump with no mid-IR emission. 
Since SCUBA observations could not be made, the SED 
could not be drawn. In contrast, the methanol maser site 
G59.78-MM2 is a much warmer clump of gas and dust. 
IRAS data was used in compiling the SED for G 59.78- 
MM2, given the close spatial coincidence between all types 
of emission. A luminosity of ~ 10^ L© and a mass of 
~ 10^ M0 were estimated for a distance of 2.1 kpc. Strong 
C^^O emission was also detected whereas the CH3CN 
lines are weak and thus only lower limits are given in ta- 
ble 5. CH3CN and C^^O fines were detected toward the 
methanol maser site at velocities falling in the methanol 
maser velocity range (16 to 27 km s"*). G59.78-MM2 is 
gravitationally bound within the mass estimate errors. 
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Table 3. CH3CN Spectral line parameters for the methanol maser sites derived from Gaussian fitting. 
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Source 


K 


Vlar 


T\ 


rms 




j' TXdv 






(km s~^) 


(K) 


(K) 


(km s"^) 


(K km s"^) 


G31.28-MM1 





110.5 


0.15 


0.01 


4.1 


0.66 




1 


110.5 


0.13 


0.01 


4.1 


0.59 
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0.24 
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0.01 
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0.14 
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0.01 


G59.78-MM2 
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21.3 
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1.8 


0.14 
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0.09 


G 173.49-MMl 
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et 
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0.10 


G 192.60-MM2 





8.2 
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0.24 




1 


7.8 


0.08 


0.02 


2.7 


0.22 




2 


7.8 


0.03 


0.02 


3.2 


0.12 




3 


7.6 


0.03 


0.02 


3.1 


0.10 



Table 4. C^^O spectral line parameters for the methanol maser sites derived from Gaussian fitting. 



Source 




Tl 


rms 




J TXdv 


A«CH3 0H 




(kms-i; 


) (K) 


(K) 


(km s-i) 


(K km s-i) 


(km s"^) 


G31.28-MM1 


108.5 


1.89 


0.17 


3.6 


7.20 


106;114 


G59.78-MM2 


22.5 


1.31 


0.07 


2.0 


2.83 


16;27 


G 173.49-MMl 


-16.7 


0.51 


0.06 


3.2 


1.71 


-15;-11 


G 188.95-MMl 


3.0 


0.66 


0.08 


3.2 


2.25 


8;12 


G 192.60-MM2 


6.2 


0.67 


0.04 


3.0 


2.15 


1;6 



3.2.3. G 173.49-h2.42 (S231, S233IR) 

G 173.49-1-2.42, also known as S231 in maser catalogues 
(e.g. Pestalozzi, Minier & Booth 2004| or as S233IR (e.g. 
Porras et al. 2000), is a methanol maser source without 
any associated visible star-forming region. The methanol 
maser site is surrounded by three Hii regions (Israel & 
Felh [WTSji . S231 (north), S233 (north-west) and S235 
(east), which are too offset to be seen in Fig. 4a. Two 
radio sources are observed east from the maser site (Fig. 
4a) and are associated with source S235 #2 in Israel & 
Felli H1978|l . No radio continuum detection has been re- 
ported toward G 173.49-h2.42 in Sridharan ct al. (511115). 
In contrast, the sub-mm continuum map is rich in emission 
sites. Four objects were detected with SCUBA: G 173.49- 
MMl, G 173.49-MM2, G 173.49-MM3 and G 173.49-MM4. 
G 173.49-MMl harbours the methanol maser, exhibits 
mid-IR emission and coincides with one of the 1.2-mm 
continuum sources observed by Beuther et al. (|2002„:) . 
G 173.49-MM2 and G 173.49-MM3 are only seen at sub- 
mm/mm wavelengths. G173.49-MM4 appears to be as- 
sociated with the far-IR source reported in the IRAS 
archives. The four mm objects are also associated with 
II2 nebulosities (Porras et al.|^00). SCUBA observations 
of G 173.49-MM3 are incomplete because it appears at 
the edge of the map. SEDs are buih for G 173.49-MMl, 
G173.49-MM2 and G 173.49-MM4, with only a single 



grey-body function fitted to the observed data of G 173.49- 
MM2. The single grey-body fit is certainly a rough esti- 
mate but it confirms that G 173.49-MM2 is the coolest 
clump. G 173.49-MMl and G 173.49-MM4 are modelled 
as massive (> 20 Mq) and luminous (> 10'^ L©) clumps 
for a distance of 1.8 kpc. The SED of G 173.49-MMl is 
not well modelled between 8.3 and 12.1 /im. This may 
be due to silicate absorption at 8.7 /im that is not taken 
into account in our modelling. CH3CN and C^^O lines 
were detected toward the methanol maser site (see tables 
3 and 4) at velocities slightly blueshifted with respect to 
the methanol maser velocity range (-15 to -11 km s~^). 
Column densities toward G 173.49-MMl derived from gas 
and dust observations agree within a factor of 2 (see ta- 
ble 5). G 173.49-MMl is gravitationally bound within the 
mass estimate errors. 

3.2.4. G 188.95+0.89 (S252, AFGL5180) 

G 188.95-1-0.89 is a methanol maser site seen in projection 
nearby Sh 2-247, a complex of Hii regions and OB stars 
(fig. 5a). Another radio source (G 188.95+0.91) to the east 
is an unknown and unresolved nebula (Kurtz et al. ll994|l . 
No radio source has been reported near the location of 
the maser site. 1.2-mm continuum emission is detected 
toward the methanol maser position. The mm map is di- 
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vided into two clumps, G 188.95-MMl and -MM2. At a 
distance of 2.2 kpc they are just resolved with the beam 
of SIMBA/SEST. G 188.95-MMl is the methanol maser 
site and exhibits emission from mid-IR to mm. G 188.95- 
MM2 is also a strong sub-mm and mm source, but does 
not coincide with any mid-IR peak emission. A third re- 
gion to the south-east appears in the MSX contour map. 
The far-IR emission observed with IRAS probably in- 
cludes emission from both G 188.95-MMl and -MM2. The 
SED of G 188.95-MMl requires a large mass (50 Mq) and 
high luminosity (~ 1.1 x lO'' Lg). The SED of G 188.95- 
MMl is not well modelled between 8.3 and 12.1 /im. This 
may be due to silicate absorption at 8.7 ^m that is not 
taken into account in our modelling. Similarly, the SED 
of G 188.95-MM2, modelled with a single grey-body func- 
tion, also requires a high mass and luminosity. Note that 
a total mass of 60 M© is given in Lada & Lada 1)2003(1 
for AFGL 5180 at a distance of 1.5 kpc. 170 Mq in total 
are derived from SIMBA observations for G 188. 95-1-0.89 
at a distance of 2.2 kpc. By taking 1.5 kpc instead of 2.2 
kpc in eq. 1, a total mass of 78 Mq would have been cal- 
culated. Our result is consistent with the Lada & Lada's 
mass estimate. CII3CN and C^^O lines were detected to- 
ward the methanol maser site (see tables 3 and 4) at ve- 
locities blueshifted with respect to the methanol maser 
velocity range (8 to 12 km s~^). Similar column densities 
are derived from gas and dust emission. The gravitational 
status of G 188.95-MMl is unclear. 

3.2.5. G 192.60-0.05 (S255IR) 

G 192.60-0.05 is a methanol maser site situated between 
two Hii regions, S255 and S257 (see Minier et al. IMITl 
and references therein), consisting of OB stars (fig. 6a). 
SIMBA observations at 1.2 mm reveal a dusty fila- 
ment elongated north-south, centred at the location of 
the methanol maser. Three mm sources are identified: 
G192.60-MM1, G192.60-MM2 and G192.60-MM3. An 
additional mm core might be present at a median po- 
sition between G 192.60-MMl and G 192.60-MM2 (Fig. 
6b), but is too poorly resolved to be successfully fit- 
ted as a compact source. G 192.60-MMl is also concident 
with a UC Hii region, a sub-mm source, and a mid-IR 
source (fig. 6b). G 192.60-MM2 is the site of the methanol 
maser. It also contains three weak and very confined ra- 
dio sources that might be interpreted as hyper-compact 
Hii regions (Minier et al. 1200 1 for explanation). It is 
also associated with mid-IR, far-IR and sub-mm emission. 
The IRAS source coincides with the maser site. G 192.60- 
MM3 is a colder source that is not clearly identified as 
a clump. G 192.60-MMl and G 192.60-MM2 are modelled 
as two spherical clumps with a high mass (> 200 Mq) 
and luminosity (> 10''' Lq) for a distance of 2.6 kpc. The 
SED of G 192.60-MMl is not well modelled between 8.3 
and 12.1 /im. This may be due to silicate absorption at 
8.7 /im. CH3CN and C^^O lines were detected toward the 
methanol maser site (see tables 3 and 4) at velocities red- 



shifted with respect to the methanol maser velocity range 
(1 to 6 km s"^). Column densities in G 192.60-MM2 de- 
rived from gas and dust observations agree within a fac- 
tor 5 (see table 5). G 192.60-MM2 is gravitationally bound 
within the mass estimate errors. 

3.3. Dust grain emissivity index ((3) 

The dust grain emissivity index (3 varies between 1.6 and 
1.9 for all sources modelled, with the exception of G 59.78- 
MM2, whose SED is poorly constrained with the absence 
of sub-millimetre data. This range of values is consistent 
with the numerical estimates of P presented by Ossenkopf 
& Henning 1)1994(1 . who found values of /3 ranging between 
1.8 and 2 for dust grains with ice mantles. The adopted 
value of Kd = 0.1 m^ kg~^ at A = 1.2 mm taken from 
Ossenkopf & Henning Iil994|l is therefore self-consistent. 



4. Discussion 

4.1. Methanol masers and MSF 

The aim of this work was to study the large-scale (~ 0.5—5 
pc) environment of methanol masers. The results demon- 
strate that each radio-quiet maser site is always associated 
with a massive (> 50 Mq) and very luminous (> 10^ Lq) 
clump of molecular gas and dust. Whether the clumps are 
very young clusters, i.e. protoclusters (< a few 10^ yr), 
forming high-mass stars or single protostellar cores is the 
subject of this discussion. 

These individual clumps are at best partially resolved 
with the relatively low angular resolution of SIMBA, 
SCUBA or MSX. An estimate of their diameter is given by 
dcoid in Fig. 2-6. It varies between 17x 10^ and 59x 10^ AU, 
i.e. ~0.1 to 0.3 pc. This implies H2 mass densities of 
6 X lO"' to 10^ Mq pc~^ in the clumps with a methanol 
maser (table 6). However, the uncertainties on the clump 
diameters is of crucial importance as the density value 
varies with d^^id- By taking the beam size of ^ 24" as 
an upper limit for each clump diameter, the mass density 
would be 5 X lO'^ to 3 x 10^ Mq pc~^ at hehocentric dis- 
tances of 1.8 to 5.6 kpc. These values agree with those 
determined for star-forming clusters such as W3:IRS5 
(Megeath et al. [TS96), Mon R2 (Carpenter et al. [TM7| . 
NGC 2024 in Orion B and the Trapezium (Lada & Lada 
12003)) . More generally, our results are consistent with the 
physical parameters describing the least evolved embed- 
ded clusters in Lada & Lada (2003) ■ Embedded clusters 
are all associated with massive (100-1000 Mq) and dense 
(riHa = 10^ — 10^ cm~'^) clumps with diameters ^ 0.5 — 1 
pc and mass density 10"^ — 10'* Mq pc^'^ within the 
molecular clouds. Interestingly, the overall mass, luminos- 
ity and linear size for each star-forming complex are not 
much larger than those measured for protoclusters forming 
stars less massive than B2 stars such as L1688 in p-Oph 
with Ltotai = 6.7 X 10^ Lq (due to HD 147889, a B2 star 
as well as the nearby Sco 0B2 association) , M ^ 500 Mq 



a. 



V. Minier et al.: Star- forming protoclusters associated with methanol masers 

^),-l°25'30" 



o 
o 
o 

C\2 



o 

o 
P 




Q -27'00 



48™00'^ 




Right Ascension (J2000) 



if" - 



G31.28-MM1 „ 


















/ ^ 






^ 

/ 


\ X 






\ <^ 






^ \ 




/ / 

p. 


\ \ 


/ 




- \ " 
' \ 




Mflas = 12.6x102 M,„,,, 


* \ 
I \ 




/To„id = 39K 






' d,„|<j = 59x103AU 


' \ 




T,„, = 110K 


' \ 
* \ 




dho, = 29 xlO^ AU 


1 \ 
1 \ 




P= 1.6 


' \ 
1 \ 



18''4B'"16" 14" 


12= 10" 


Right Ascension 


(J3000) 


G31.28-MM2 








X ^.v ^ 

/ Mgas = 3.4x102|yl3'„,,, v 
T,„„ = 16K 

d,„id = ioaxio3AU 


\ 

\o 

\o - 


T,„,= 145K 

df,o^=6x102AU ' 


\ 

\ 

\ 

\ 


43'= 1.6 i_ 


\ 

\ 

. . \ 



Frequency (GHz) 



io' 10' 

Frequency (GHz) 



10' 



d. 





150 

Velocity [km/s] 



100 no 

Velocity [km/s] 



Fig. 2. G 31.28+0.06. a. Optical image (grey scale) of the neighbourhood of G 31.28+0.06. The grey contours represent 

the 1.2-mm contimmm emission imaged by SIMBA (5, 10, 20, 50 and 90% of peak flux). The cross represents the 
position of the methanol masers, the squares are the radio continuum sites, the diamond is a recombination line site 
and the large circle shows the position of the IRAS source. The methanol maser source is not coincident with any 
visible object or optically visible Hii region but it is associated with a strong mm source, b. Closc-up of the mm 
sources 1 and 2. 850-/im SCUBA (grey scale), 450- /xm SCUBA (grey contours and 10, 50 and 90% of peak flux) and 
21-/um MSX (white contours and 20 and 60% of peak flux) images. The mm sources 1 and 2 are also seen with SCUBA 
and MSX. c. SEDs of the mm sources 1 and 2. The IRAS icon shows the IRAS 100-/um flux level, d. CH3CN and 
C^^O line spectra. 
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Fig. 3. G 59.78+0.06. a. Optical image (grey scale) of the neighbourhood of G 59.78+0.06. The grey contours represent 

the 1.2-mm contimmm emission imaged by SIMBA (10, 20, 50 and 90% of peak flux). The cross represents the position 
of the methanol masers, the squares are the radio continuum sites and the large circle shows the position of the IRAS 
source. The methanol maser source is not coincident with any visible object or optically visible Hll region but it is 
associated with a strong mm source, b. Close-up of the mm sources 1 and 2. 1.2-mm SIMBA (grey scale and contours; 
10, 20, 50 and 90% of peak flux) and 21-/xm MSX (white contours; 20 and 60% of peak flux) images. The mm source 
2 is also seen with MSX. c. SED of the mm source 2. d. CH3CN and C^^O line spectra. 
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Fig. 4. G 173.49+2.42 (S231). a. Optical image (grey scale) of the neighbourhood of G 173.49+2.42. The grey contours 

represent the 850-/im continuum emission imaged by SCUBA (10, 20, 50 and 90% of peak flux). Note that the SCUBA 
map does not entirely cover the optical map. The cross represents the position of the methanol masers, the squares 
are the radio continuum sites and the large circle shows the position of the IRAS source. The methanol maser source 
is not coincident with any visible object or optically visible Hii region but it is associated with a strong mm source, 
b. Close-up of the mm sources 1, 2, 3 and 4. 850- /um SCUBA (grey scale), 450- /xm SCUBA (grey contours; 10, 20, 50 
and 90% of peak flux) and 21-/im MSX (white contours; 20 and 60% of peak flux) images. The mm sources are also 
seen with SCUBA, and MMl and MM4 are seen with MSX. c. SEDs of the mm sources 1, 2 and 4. The IRAS and 
MSX icons show the IRAS 100-/xm flux level and MSX 20-/xm image noise rms level, d. C^^O line spectrum. 
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Fig. 5. G 188.95+0.89 (AFGL5180). a. Optical image (grey scale) of the neighbourhood of G 188.95+0.89. The grey 

contours represent the 1.2-nim continuum emission imaged by SIMBA (10, 20, 30, 50. 70 and 90% of peak flux). 
The cross represents the position of the methanol masers, the squares are the radio continuum sites and the large 
circle shows the position of the IRAS source. The methanol maser source is not coincident with any visible object 
or optically visible Hii region but it is associated with a strong mm source, b. Close-up of the mm sources 1 and 2. 
850-/im SCUBA (grey scale), 450-/xm SCUBA (grey contours; 30, 40, 50, 60, 70 80 and 90% of peak flux) and 21-/im 
MSX (white contours; 20 and 60% of peak flux) images. The mm sources 1 and 2 are also seen with SCUBA and MSX. 
c. SEDs of mm sources 1 and 2. The IRAS icon shows the IRAS lOO-yum flux level, d. CH3CN and C^^O hue spectra. 
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Fig. 6. G 192.60-0.05 (S255IR). a. Optical image (grey scale) of the neighbourhood of G 192.60-0.05. The grey contours 

represent the 1.2-mm contimmm emission imaged by SIMBA (10, 20, 50 and 90% of peak flux). The cross represents 
the position of the methanol masers, the squares are the radio continuum sites, the diamond is a recombination line 
site and the large circle shows the position of the IRAS source. The methanol maser source is not coincident with 
any visible object or optically visible Hii region but it is associated with a strong mm source, b. Close-up of the mm 
sources 1 and 2. 850-/Ltm SCUBA (grey scale), 450- /um SCUBA (grey contours; 25, 45, 65 and 85% of peak flux) and 
21-/im MSX (white contours; 5, 45 and 85% of peak flux) images. The mm sources 1 and 2 are also seen with SCUBA 
and MSX. c. SEDs of mm sources 1 and 2. The IRAS icon shows the IRAS 100-/im flux level, d. CH3CN and C^^O 
line spectra. 
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Table 5. Physical properties of the methanol maser sites derived from the SED and spectral line analysis. Nn^ and 
riHa are derived from Mgas in the continuum case by using an angular size of 24". The values of Nn^ derived using 
dcoid are given in parentheses. Mgas in the spectral line case is derived from Nn^ by using an angular size of 34". 
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G 31.28-MMl 
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12.6 


17.5 (90.3) 


13.1 


10.3 


34 


21.6 


15.1 


8.0 


7.4 


G59.78-MM2 


53 


1.0 


9.8 (36.5) 


19.5 


1.1 


> 18 


> 0.8 


> 3.8 


> 5.4 


0.9 


G 173.49-MMl 


49 


1.2 


15.5 (31.8) 


36.0 


4.9 


69 


1.0 


6.5 


10.8 


1.8 


G 188.95-MMl 


42 


0.5 


4.5 (41.0) 


8.6 


1.1 


58 


1.6 


7.4 


10.0 


2.3 


G 192.60-MM2 


42 


3.2 


22.3 (59.5) 


37.3 


5.1 


55 


1.5 


5.3 


6.3 


2.4 



and a diameter ~ 1 pc (Greene & Young EHU Wilking & 

Ladannnni). 

The visual extinction of the clumps with a methanol 
maser can be estimated with Ay = N-h_.^/Q.QA x 10^^ mag 
(Frerking et al. 11982(1 . It varies from 40 to 240 mag, im- 
plying a high degree of embeddedness and presumably a 
very young age (<C 10^ yr). In comparison, the visual ex- 
tinction in the p-Oph central C^^O cloud in L1688 is in 
the range 50-100 mag (Wilking & LadaHlHSI) with a col- 
umn density ~ 10^^ cm~^ and a diameter ~ 0.5 pc. Three 
mm cores (Oph-C, Oph-CS and Oph-E) were mapped by 
Motte et al. p998|l in the central cloud and have column 
densities ranging from ^ 7x 10^^ to 10^'^ cm~^, which com- 
pare to those in the methanol maser clumps (4.5 x 10^^ to 
2.2 X 10^3 cm-2). 

In summary, the methanol maser sites exhibit the 
properties of young, deeply embedded clusters of YSOs. 
Only luminosities in excess of 10"* L0 (^^ZAMS BO. 5, 
Panagia I1973|l for clumps less than ~ 1 pc in diameter 
clearly suggest the presence of high-mass (> 8 M©) YSOs. 

The detection of CH3CN and C^^O toward each 
methanol maser site indicates that they are within HMCs 
with gas density > 10^ cm"^. CH3CN emission usually 
arises from the 0.1-pc inner part of the hot core although 
ground state CII3CN emission could trace a cooler region 
in HMCs than those probed with vibrationally excited 
CH3CN lines (Olmi et al. 2003). This might explain the 
relatively low rotational temperatures (18-69 K) in table 5. 
The C^*0 antenna temperatures (table 4) are much larger 
than those measured for CII3CN (table 3). A possible ex- 
planation is that C^^O emission arises from a larger area 
than CH3CN emission. The CH3CN rotational tempera- 
tures agree well with the Tcoid values (table 5), which sug- 
gest that CH3CN(5if — 4^) and dust emission arises from 
the same area. The fitted source solid angles {^coid and 
^hot used to estimate dcoid and dhot) in the dust emission 
models are all smaller than the OSO and Mopra telescope 
beam sizes, which then probably overestimate the true 
size of the emitting regions. In consequence, the densities 
(~ 10^ cm-3) derived from the C^^O data (table 5) with 
the telescope beam sizes as the source angular sizes are 
lower limits of the CII3CN clump densities. In addition to 
CII3CN emission, NII3 and CII3OII quasi-thermal lines 



were also detected toward G 59.78-f0.06, G 173.49-h2.42, 
G 188.95-^0.89 and G 192.60-0.05 (Minier & Booth I^HH^ . 
The presence of the typical HMC species (CII3CN, NH3 
and CII3OH) suggest again a very young age (< a few 10^ 
yr) for the methanol maser sites (e.g. Rodgers & Charnley 
(2001). 

In conclusion, methanol masers trace very luminous, 
deeply embedded, massive, molecular clumps in the 
five star-forming regions considered here. In a few of 
the clumps (G 31.28-MMl, G 192.60-MM2 and G 59.78- 
MM2), ionised sources were also detected in the radio. 
This implies that high-mass stars are currently forming 
within these clumps and that the deeply embedded star- 
forming clumps are protoclusters of high-mass YSOs. 

Yet, whether methanol masers are themselves di- 
rectly associated with a massive YSO remains unclear. 
In G31.28-h0.06, the maser is offset by 3.5 x lO"* AU 
from the UC Hli region (Walsh et al. 1998). The maser 
site in G 59.78-1-0.06 appears to coincide with a very 
weak radio source within 2" (Sridharan et al. I2()()2(l . 
In contrast, no radio continuum emission has been de- 
tected in G173.49-h2.42 and in G 188. 95-f 0.89. Finally, 
the methanol maser in G 192.60-0.05 coincides with a 
4.4-mJy radio source. Assuming that a maser site and 
a radio continuum source coincident in the maps are 
indeed co-spatial, the methanol masers in G 59.78-1-0.06 
and G 192.60-0.05 are associated with the precursors of 
a UG Hii region. In G 31.28-1-0.06, G 173.49-^2.42 and 
G 188.95-1-0.89, further high angular resolution and high 
sensitivity mm continuum observations are required to 
identify the nature of the maser sources. 

4.2. Multiple star-forming (proto)clusters 

A secondary aspect of this work is the detection of objects 
only seen in mm continuum (hereafter "mm-only" clumps) 
and the variation of the SED profiles with the evolutionary 
status of the YSOs. 

4.2.1. Cold cores 

G59.78-MM1 and G 173.49-MM2 are typical mm-only 
clumps. G 31.28-MM3, G 173.49-MM3 and G 192.60-MM3 
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Fig. 7. Luminosity versus mass and colour diagrams. 
Each number represents an object as follows G 31.28- 
MM1(#1), G31.28-MM2(#2), G59.78-MM2(#3), 
G173.49-MM1(#4), G 173.49-MM2(#5), G 173.49- 
MM4(#6), G188.95-MM1(#7), G 188.95-MM2(#8), 
G192.60-MM1(#9), G 192.60-MM2(#10). The dashed 
line in the top plot represents Ltotai^Mg^^ . The dashed 
lines in the colour diagrams represent a black body 
function with the temperature ranges indicated on the 
plots. The dotted line in the middle plot represents the 
massive YSO colour criterion from Lumsden et al. (EE 



are other possible mm-only clumps, making in total 
five mm-only clumps. G 188.95-MM2 probably emits in 
mid-IR, but its flux was not measured due to con- 
fusion with the main mid-IR source, G 188.95-MMl. 
These mm-only clumps are cold (~ 20 K for G 173.49- 
MM2). The modelled luminosity in G173.49-MM2 is 
only 300 Lq. Assuming Tcou = 20 K for G31.28-MM3, 
G59.78-MM1, G 173.49-MM3, G 192.60-MM3, masses of 
~ 25 — 178 M0 are derived for the five mm-only clumps. 
Excluding G 31.28-MM3 at 5.6 kpc, masses range 19 — 
103 Mq for G 59.78-MMl, G 173.49-MM2, G 173.49-MM3 
and G 192.60-MM3. The modelled diameter of G 173.49- 
MM2 is only 21 x 10^ AU, i.e. 0.1 pc. Using an angular 
size of 24" for each object, column densities ^ 2.5 — 7.4 x 
10^^ cm~^ and visual extinction ~ 26 — 78 mag are esti- 
mated. The physical parameters of all detected mm-only 
clumps are summarised in table 6. 

The cold clumps might be prestellar clumps or 
protostellar clumps in the process of forming stars. 
Unfortunately, such types of clumps have marginally been 
detected and it is difficult to compare the mm-only clump 
properties with those of statistical samples of massive 
prestellar and early protostellar clumps. A good protostar 
candidate in an early accretion phase is IRAS 23385-1-6053 
(Fontani et al. I2003p . which is a cold (~ 40 K) and small 
(~ 0.03 pc) molecular core embedded in a colder (~ 15 K) 
and larger (< 0.4 pc) halo. The SED, luminosity (150 L©) 
and temperature (15 K) of IRAS 23385-1-6053 resemble 
those of G 173.49-MM2, when the IRAS points are ex- 
cluded (Fig. 13 in Fontani et al. I2003|) . But as stated by 
the authors, it represents a lower limit. The luminosity of 
G 173.49-MM2 could also be underestimated without the 
knowledge of the far-IR flux. 



4.2.2. SED profiles 

In addition to the mm-only clumps, three other types of 
SED profile are identified. They might represent different 
evolutionary phases or result from geometrical effects that 
modify our perception of the YSO environments in func- 
tion of the linc-of-sight. G31.28-MM2 is the only example 
of a cold (16 K) envelope with a nearly hot (115 K) core. 
The envelope continuum emission and the core continuum 
emission peak at relatively low frequencies. G 59.78-MM2, 
G 173.49-MM4, G 188.95-MMl and G 192.60-MM2 exhibit 
SEDs with a "warm" (~ 45 K) component peaking in the 
far-IR at 60 — 100 ^vci and a hot (150 — 250 K) component 
peaking in the mid-IR at 15 — 30 /im. They are all very 
luminous (> 10"* Lq). Finally, G31.28-MM1, G173.49- 
MMl, G192.60-MM1 and perhaps G 188.95-MM2 see their 
SED profile dominated by a strong far-IR component. The 
core component of the SED has shifted to longer wave- 
lengths (~ 50 ^m) and cooled down to ^ 110 K. The 
warm (~ 40 K) envelope emission has strengthened and 
still peaks around 100 /im. The temperatures of all the 
"warm" components are in good agreement with those ex- 
pected given the Imirinosities, sizes and (3 (e.g. see eq. 11 
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in Garay & Lizano [T999|l . The hot component most likely 
probe the central part {dhot = 100 to 3000 AU in Fig. 2-6) 
of the core around the luminous protostellar object. 

To compare the evolutionary status of each ob- 
ject, arguments are presented below and in Fig. 7. 
Only the clumps for which the SED can be built 
are considered. G173.49-MM2 and G31.28-MM2 present 
Lsubmml Ltotai > 10^^ while the other clumps have 
Lsubmm/Ltotai < 3 X 10"^, where Lsubmm is the luminos- 
ity integrated between 1.2 mm and 450 fim. Similar results 
are noticed when considering Ltotai 1^1^ gas ratios, which are 
also independent of distances. G 173.49-MM2 and G 31.28- 
MM2 have lower Ltotai/Mgas ratios (3-10 Lq Mq^) than 
those measured for the other clumps (> 80 L© Mq^). 
This could be explained if G173.49-MM2 and G 31.28- 
MM2 are younger than the other objects associated with 
maser sites and UC Hii regions, hence less luminous. This 
is also seen in Fig. 7-top where G 173.49-MM2 (#5 in 
Fig. 7) and G31.28-MM2 (#2 in Fig. 7) lie below the 
other sources whose luminosities follow Ltotaio^Mg-J^ . The 
objects G 59.78-MM2, G 173.49-MM4, G 188.95-MMl and 
G 192.60-MM2 are all relatively hot inside, and strongly 
emit in the mid-IR with respect to the sub-mm emis- 
sion (Fig. 7 bottom). They are all good HMC candidates. 
The mid-IR component is very much absorbed in the 
FIR component in the SEDs of G31.28-MM1, G 173.49- 
MMl, G188.95-MM2 and G192.60-MM1. G31.28-MM1 
and G192.60-MM1 contain Hii regions. The cases of 
G 173.49-MMl and G 188.95-MM2 are unclear and they 
might belong to the HMC or UC Hii class of objects. 

Alternatively, the least massive clumps near the 
methanol maser site could be low- or intermediate-mass 
star-forming protoclusters. Their physical characteristic 
resemble those of the p-Oph mm cores (Motte et al. 11998(1 . 
Indeed G 173.49-MM4 is coincident with S233IR NE, a 
cluster of low-mass YSOs seen in near-IR (Porras et al. 
Onnn|l . The derived luminosity (6 x 10^ L©) and mass (20 
Mq) are typical of the mm cores in p-Oph. In Fig. 7- 
middle, G 173.49-MM4 appears slightly isolated from the 
group of clearly identified massive star-forming clumps 
for which Fu > Fg and F2i/F% > 2 (e.g. Lumsden et 
al. l20U2fl . G 173.49-1-2.42 is therefore a good example of a 
star-forming complex in which various mass types of YSO 
clusters co-exist. 

High resolution and high sensitivity observations in 
mid-IR, mm and cm are needed to differentiate between 
these two scenarii by searching for multiple mm core sys- 
tem within each mm clump as seen in p-Oph and weak 
free-free emission around high-mass protostars. 

5. Conclusions 

Five methanol maser sites have been studied from mid- 
IR to mm wavelengths. Each radio-quiet maser site is al- 
ways associated with a massive (> 50 M©), very lumi- 
nous (> 10^ Lq), deeply embedded (> 40 mag) molec- 
ular clump. They are all underluminous stellar objects 
[Ltotaio^Mg-Jg ). The properties of the maser sites derived 



from observations and modelling suggest that they are em- 
bedded HMCs in which high-mass stars are in the process 
of forming. All these elements demonstrate that the radio- 
quiet methanol masers trace high-mass star-forming pro- 
toclusters in earlier phases than those detected in the H ii 
region surveys. 

In addition, colder gas clumps seen only at mm wave- 
lengths are also found near the methanol maser sites. 
These might represent an even earlier phase of massive 
star formation. These results suggest an evolutionary se- 
quence for massive star formation from a cold clump 
only seen at mm wavelengths, evolving to a HMC with 
a two-component SED with peaks at far-IR and mid-IR 
wavelengths, and ending with (ultra-compact) H ii regions 
very bright at far-IR wavelengths. Alternatively, the cold 
clumps might be clusters of less massive (< 8 M©) YSOs, 
in formation near the high-mass star-forming protoclus- 
ters. 

Finally, the values of the dust grain emissivity index 
(/?) range between 1.6 and 1.9, which agree with modelling 
work (Ossenkopf & Henning .1994;) . 
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